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I Abstract 

ASI Inspired by the recent LHC Higgs data and null search results of supersymmetry (SUSY), we 
eS scan the parameter space of the Minimal Supersymmetric Standard Model (MSSM) with relatively 
TT ; heavy sparticles (1-3 TeV). Then in the parameter space allowed by current collider experiments 
c» i; and dark matter detections, we calculate the complete one-loop SUSY QCD corrections to pp — bbh 
oO at the LHC with ys = 14 TeV and obtain the following observations: (i) For the large values of 
(Com 

Í* à tan f and low values of m4, the SUSY QCD effects can be quite large, which, however, have been 
Cy t 

‘AL: excluded by the latest results of LHC search for H/A — 7*77; (ii) For modest values of tan 3 and 
i e. mA which so far survived all experimental constraints, the SUSY QCD corrections can maximally 


reach about — 996. 


PACS numbers: 


I. INTRODUCTION 


Very recently the ATLAS and CMS collaborations have independently reported the ob- 
servation of a Higgs-like resonance with a mass about 125 GeV [I]. At the same time, the 
CDF and D0 collaborations have also updated their combined results for the Higgs searches 
in bb channel, which support the LHC observation B. Since in the Minimal Supersymmetric 
Standard Model (MSSM) a SM-like Higgs boson is predicted with a mass below 130 GeV, 
the observation of such a 125 GeV Higgs boson supports SUSY, albeit quite restrictive on 
the parameter space of SUSY |3|. 

Meanwhile, the direct searches for SUSY particles (sparticles) have been performed at the 
LHC. Based on about 5 fb^! luminosity, the ATLAS and CMS collaborations have reported 
null results and obtained some bounds on the sparticle masses, which is about 1 TeV for the 
gluino and first generation of squarks |4], 330 GeV for the electroweak gauginos, 180 GeV 
for the sleptons li 465 GeV for the stops and 480 GeV for the sbottoms [6]. These bounds 
indicate that SUSY may be heavier than expected and the sparticles may be significantly 
heavier than the electroweak scale |7, 

In case that the sparticles are heavy and beyond the LHC scope of direct production, 
search for the indirect SUSY effects via loop corrections will be of great importance. Since 
the loop effects of heavy sparticles are usually small, we should look for some processes 
in which the heavy sparticles have residual loop effects. One type of such processes are 
Higgs productions at the LHC, such as the production of tH~ and Abb, in which the heavy 
sparticles have sizable residual loop effects for a small value of m, and a large value of tan 6 
i 10] (when m4 getting large, such effects will vanish). The reason for these residual loop 
effects is that the couplings in the loops are proportional to some SUSY mass parameters 
and can be enhanced by the large values of tan f. 

In this note we focus on the production of Abb at the LHC and calculate the complete one- 
loop SUSY QCD corrections to this process. As an important Higgs production channel for 
the MSSM, this production has been studied in the literature , where the residual SUSY 
QCD effects are found to be large (reach -40% for tan 3 = 30). We revisit this production 
for the following reasons: Firstly, in the literature the SUSY QCD corrections to this process 
are calculated only partially (only the corrections to the hbb vertex have been considered). 


The complete one-loop corrections involve pentagon Feynman diagrams, whose calculations 


are rather complicated and have not been performed. Secondly, the CMS collaboration has 
recently measured this channel and given constraints on the plane of tan 9 versus mA . 
Since the residual SUSY QCD effects in this production is sensitive to the values of tan 6 
and m4, we should update the calculations by considering such new constraints. Moreover, 
other experimental constraints, such as the dark matter direct detection limits and the SM- 
like Higgs boson mass around 125 GeV, are also rather restrictive and should be considered. 
In this work, we consdier all current experimental constraints to scan the MSSM parameter 
space and then in the allowed parameter space we calculate the process pp — bbh with the 
complete one-loop SUSY QCD corrections. 

The paper is organized as follows. In Sec. II. we describe the calculations for the process 
pp — bbh. In Sec.III we show numerical results. Finally, we draw the conclusions in Sec. 


IV. 


Il. THE DESCRIPTION OF CALCULATIONS 


In the MSSM the lighter CP-even Higgs mass (mj) is smaller than Mz at tree level but 
receives large corrections at the loop level. The leading part of the corrections is from the 


stop sector and can be expressed as |13 


5 3m4 Mz m; X? xX? 
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where X, = A,—p cot B is the stop mixing parameter. We see that a large stop mass or a large 
stop mixing parameter is needed to increase mp, to 125 GeV. In our calculations we consider 
the collider constraints on the MSSM Higgs sector, using the packages FeynHiggs2.8.6 
and HiggsBounds-3.8.0 to calculate the observables in the Higgs sector and require them 
to satisfy the constraints from the LEP, Tevatron and LHC. 

The SUSY QCD corrections to hbb production at the LHC involve the sbottoms and 


gluino in the loops. The sbottom mass matrix takes the form |16 


where 


1 1., 
m; = mo +m -m (5 ^3 sin” Oy) cos(26) , (3) 
1 
y 0e Magog raz ein! Pure (aia » 
Xy = Ay — ptanf , (5) 


where ms and m are respectively the soft-breaking mass parameters for the left-handed 


squark doublet Q and the right-handed down squark D, A, is the sbottom soft-breaking 
trilinear coupling and p is the SUSY-preserving bilinear coupling of the two Higgs doublets 
in the superpotential. This mass matrix can be diagonalized by a unitary transformation 


which rotates the weak eigenstates br, r to the mass eigenstates bi 2, 


b, cos; sin6; b L 


b; —sin 6; cos 6; bp 


with the sbottom masses m; , and the mixing angle 6; determined by 


1 2 
DR. = 5 m + my T (mj, - m.) + m , (7) 
2myX 
Mir Mir 


We produce the one-loop amplitudes with FeynArts-3.5 and use the FormCalc-6.1 
to simplify them and express the loop functions. The numerical calculations are performed 
by using LoopTools-2.2 |19]. In Fig.1 we display the representative pentagon Feynman 
diagrams for the SUSY QCD corrections in the subprocesses gg — bbh. Due to no massless 
particles in the loop, all the Feynman diagrams with the gluino and sbottoms in the loops 
are infrared (IR) finite. 

We take the definitions of the scalar and tensor two-, three-, four- and five-point inte- 
gral functions presented in Ref.|20]. For the calculation of the pentagon diagrams, we use 
Passarino-Veltman method|21] to reduce the N-point(N < 5) tensor functions to scalar in- 
tegrals. Our programs have been used to study the SUSY-QCD corrections to the process 
pp — ttZ? at the LHC and have been checked with Mn therein. In order to further 
validate the calculation of the pentagon diagrams, we used our programs to calculate the 
NLO QCD corrections to pp > tth in the SM at the LHC and compared with the results in 
net od 


. As shown in Table I, our results argee with those in Ref. very well. In order to 
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(11) 
FIG. 1: The pentagon diagrams for SUSY QCD corrections to the subprocess of gg — bbh at the 


LHC. 


preserve supersymmetry, we adopt the constrained differential renormalization (CDR) 
to regulate the ultraviolet divergence (UV) in the self-energy and vertex corrections, which 


is equivalent to the dimensional reduction method at one-loop level . 


TABLE I: The comparison of our numerical results for the process pp — tth in the SM at the 
LHC with those in nld where the LO and NLO QCD corrected cross sections for different 
Higgs mass are listed with the relevant parameters and the PDFs being the same as in Ref.|24], 


i.e., ui = (2m, + m4)/2, m, = 174 GeV and the MRST PDFs. 


INC vol) Bi m e 
TA 
(12 


xa xs 
373.6(2) | 4524(11) 373.4(3) 452.3(12) 
251.3(4) | 305.5(7) 251.6(2) 305.6(8) 


In our calculations, we assume a common SUSY mass Msysy defined by Msysy = Ms = 


Mg = Mg = Mg = A, = Ay = u. We fix slepton mass parameters M; = Mg = A, = 3 TeV 


and scan the following MSSM parameter regions: 
5 <tan6<60, 90GeV < M4 € 350 GeV, 1 TeV € Msusy €3 TeV (9) 


In our scan we consider the following constraints on the parameter space: (1) We require that 
the bounds for Higgs bosons from LEP, Tevatron and LHC are satisfied and the mass of light 
CP-even Higgs is in the region of 123 GeV « mp < 129 GeV; (ii) For the constraints from 
flavor physics and electroweak precision tata, we checked by using the package susy.flavor 
v2.0 " that they are safely satisfied because we assume relatively heavy sparticles. (iii) We 


consider the dark matter constraints from the WMAP relic density and the direct detection 


results by using the package MicrOmega v2.4 : 


Ill. NUMERICAL RESULTS 


Since the b-quark Yukawa coupling may receive large radiative corrections in the MSSM, 
we use the running b-quark mass (mP®) and use the method induced in to absorb the 
MSSM corrections into the effective b-quark Yukawa couplings. But for the b-quark in the 
final state, we take the pole mass to assure the correct on-shell behavior. 


In our numerical calculations, we take the input parameters of the SM as 


m, = 172 GeV, m¥5(m™$) = 419 GeV, mz = 91.1876 GeV, a(mz) = 1/127.918 


Here mI? (mM) is the QCD-MS bottom-quark mass, which is related to mp” as 


mPF = qfs[ — 2 o (73 — 3n] (10) 
=m m nales — 3n 
b 3r 144r? : 


where n; is the number of the active quark flavors. For the strong coupling constant o;(11), 
we take its 2-loop evolution with QCD parameter A"/-? = 226 MeV and get a,(mz) = 0.118. 
We use CTEQ6L1 and CTEQ6M parton distribution functions (PDF) for the SM tree 
level and SUSY QCD one-loop level computations, respectively. The renormalization scale 
ur and factorization scale up are chosen to be ur = ur = mz. We numerically checked 
that all the UV divergence in the loop corrections canceled. 

In Fig.2 we project the survived samples satisfying all the experimental constraints on 
the planes of tan 9 and Mp versus m4 and also present the excluded regions. It can be seen 


that a large part of the parameter space (the light blue region) has been ruled out by the 
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FIG. 2: The parameters space satisfying constraints (i)-(iii), projected on the planes of tan 3 versus 
ma. The blue region is excluded by the LHC data, in particular, by the search of new particles 
decaying into 7T* 7^ and the measurement of B, > utp. The yellow region is excluded by the 
non-observation of Higgs boson at LEP2. The red dots represent the samples survied all the 


constraints. 


7 TeV LHC data, in particular, by the search of new particles decaying into 7*7 ^ and the 
measurement of B, > u^. For the small tan 8 and low and moderate M4 region, it has 
been excluded by the non-observation of Higgs boson at LEP2. We also note that, with the 
very recently released 7+8 TeV LHC results of H/A > r* r^ based on £ = 17 fb [33], the 
excluded lower limit on the plane of tan 9 — M4 has been further pushed down and overlaps 
with the one of LEP2 in the low M4 case. Since we require the Higgs boson mass to be 
in the range of 123-129 GeV indicated by the LHC data(126+3 GeV), the parameter space 
that can correctly produce the Higgs boson mass is highly constrained and situated in a 
region with a modest M4 (M4 Z 200 GeV) and a small tan 8 (6 X; tan 8 < 12). However, 
for other parts of the parameters space, they produce a too heavy Higgs boson (mj, > 129 
GeV when tan B > 12) or a too light Higgs boson (mp < 123 GeV when tan 8 < 6). 

In Fig.3 we present seperately the pentagon diagram contribution (lower panel) and 
the total SUSY-QCD contribution (upper panel) for the surviving samples. In order to 
show the influence of the recent LHC data, we present the complete one-loop SUSY QCD 
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FIG. 3: The upper panel shows the complete one-loop SUSY QCD corrections(da/o) versus MA 
and tan( for the samples satisfying(red triangles) or not satisfying(green dots) the LHC con- 
straints. The lower panel shows the SUSY-QCD corrections from the contributions of pentagon 


diagrams(d0,/o) for the red triangle samples. 


corrections for the samples satisfying or not satisfying the LHC constraints. We can see 
that the complete SUSY QCD corrections will be significant for the samples which have 
a large tan and a low value of m4. This can be understood by the contribution to the 


effective b-quark Yukawa coupling after integration of the heavy sparticles, which is Ypa = 
Em S7 (Me tan B + cota] [8]. Since we assume Msusy = M; = u, the Yukawa 
coupling will be independent of the sparticle masses and be greatly enhanced by a large tan f. 
However, it should be noted that these samples will lead to the excess of the production rate 
of pp > H/A — r*7- and thus have been excluded by the current measurements. With the 
increase of ma and the decrease of tan £, the corrections drop rapidly and approach zero in 
the decoupling limit. The main reason is that Y, can be heavily reduced by the cancellation 


between tan and cosa, which have a relation as cota ~ — tan 8 — 2m? tan B cos? B /m?, 


for a large m4. From the lower panel of Fig.3, we see that contributions of those pentagon 
diagrams are small and maximally reach about 0.896. 'This is because in those pentagon 
diagrams the Higgs boson only couples to sbottoms while the Abb vertex does not appear. 
So the large residual loop effects in the hbb vertex are absent in the pentagon diagrams. For 
the samples which survived all the constraints, the complete SUSY QCD corrections can 
only reach about —9% at the LHC with ys = 14 TeV. Detecting such a size of SUSY QCD 


effects may be challenging in the future measurement of the process pp — bbh |34]. 


IV. CONCLUSION 


In this work, we calculated the complete one-loop SUSY QCD corrections to the process 
pp — bbh at the LHC with \/s = 14 TeV. We found that the large SUSY QCD corrections in 
the non-decoupling regime with a large tan 6 and a low my has been excluded by the latest 
results of LHC non-standard Higgs searches. For the survived decoupling regime which have 


modest values of tan 8 and m4, the SUSY QCD corrections can maximally reach —9%. 
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